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ABSTRACT 

We use numerical simulations to study the kinematic structure of remnants formed from mergers of 
equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include 
the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive 
black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show 
significant rotation compared to remnants of dissipationless mergers. The increased rotation speed 
of dissipational remnants owes its origin to star formation that occurs in the central regions during 
the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis 
rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more 
isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky 
than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A 
large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all 
of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers 
can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/cr)*, 
kinematic misalignments, and isophotal shapes. The dissipationless remnants are a poor match 
to this data. We also investigate the properties of merger remnants as a function of initial disk gas 
fraction, orbital angular momentum, and the mass of the progenitor galaxies. Our results support 
the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor 
disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that 
are slowly rotating and uniformly boxy. 

Subject headings: galaxies: ellipticals — evolution — formation — interactions — methods:N-body 
simulations 



1. INTRODUCTION 

The observed absorption-line spectra and red colors of 
elliptical galaxies suggest that their stars were formed 
at high redshift (z > 1) and that very little star forma- 
tion has occurred in t hem since then. According to the 
me rger hypothesis" l|Toomre fc Toomreill97l iToomrel 
119771) . these red elliptical galaxies are produced by the 
collision and merger of spiral galaxies, and hence the pro- 
genitors of present day ellipticals may be high-redshift 
spirals. While relatively little is known about disk galax- 
ies at high redshift, it is likely that these disks were more 
concentrated and gas-rich than their low-redshift coun- 
terparts. Indee d, preliminary observational evidence 
ijErb et al.l2006[) indicates that galaxies at redshift z w 2 
do have large gas fractions / gas ~ 0.5, with some ap- 
proaching / gas ~ 0.8 — 0.9. Thus, any attempt to under- 
stand the formation, properties and scaling relations of 
the present day population of elliptical galaxies, within 
the context of the "merger hypothesis" , must consider 
gas-rich mergers. 

Requiring that the disk galaxy progenitors contain a 
significant fraction of gas is nothing new to the "merger 
hypothesis". One of the main objections to this mech- 
anism of prod ucing ellipticals, argued by, for instance, 
lOstrikerl Jl980), is that the centers of ellipticals are more 
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concentrated than local spirals. Cast in terms of phase- 
space density, this objection states that the high central 
phase-space density of ellipticals cannot be produced by 
the merger of low phase-space spirals because, according 
to Liouville's Theorem, phase -space density i s conserved 
during a collisionless process l|Carlberdll986|) . However, 
this argument breaks down when the merger constituents 
contain gas, which can r adiate energ y, and hence increase 
the phase space density (|Lakell989j) . An estimate of how 
much gas is required to match the central densi ties of 
ellipticals was provided bv lHernauist et aD 1^993), who 
used N-body simulations and analytic arguments to sug- 
gest that mergers of spiral galaxies containing > 30% gas 
would be sufficient to account for the high phase space 
densities of ellipticals. 

Within the context of the hierarchical theory of struc- 
ture formation, gas-rich major mergers may play a 
much larger role than just resolving the central phase- 
space densities of elliptical galaxies. Previously, we 
have described a "cosmic cycle" of galaxy formation 
and evolution in w hich gas-rich mergers drive the evo- 
lution of quasars (|Hopkins et alJ l2006al). induce the 
grow th of supermassive black holes |Di^fattej3_e_t_al ' 



2005J), and produce red ellipt ical galaxies l|Springei"et*ld 



2005al IHopkins et alJl2006cD that obey many of th e ob- 
served scaling relations l|Robertson et alJ l2006albl) . A 

schematic view of thi s picture is presented in Fig. 1 of 
IHopkins et al.l ( 2006a). However, the success of this sce- 
nario also must be gauged by its ability to produce rem- 
nants that have kinematic and morphological properties 
characteristic of observed ellipticals. It is this question 
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that we address in the current paper. 

Observations indicate that gal axy spheroids can 
be classified into two groups llDavies et all 1983; 
Bender et alJ Il989t iBenderl [1983 iFaber et all 11997 : 
KOTm^nTivfcBenderl 119961 and references therein). 
Large, luminous spheroids have hot gaseous halos, box- 
shaped isophotes, surface-brightness profiles with flat 
"cores," show very little rotation and are almost uni- 
formly classified as ellipticals. Less luminous spheroids 
tend to have little, if any, hot gas, disk-shaped isophotes, 
power-law surface-brightness profiles, and exhibit rota- 
tion that is along the photometric major axis. The lat- 
ter group encompasses many low-luminosity ellipticals, 
bulges and SOs. 

The last of these properties, the alignment and rota- 
tion of spheroidal galaxies, may be the result of a more 
fundamental dichotomy among elliptical galaxies; that 
is, the isotropy of their velocity distributions. Because 
spheroids can be flattened for (at least) two reasons, ro- 
tation and velocity anisotropy, the lack of rotation in 
large ellipticals suggests that these systems have signif- 
icant velocity anisotropy, while many low-luminosity el- 
lipticals and bulges are consistent with being isotropic 
systems flattened by their observed rotation. 

One viewpoint is that these two classes of elliptical 
galaxies exhibit different properties because they are 
formed yia different mechanism s. Along this line of rea- 
soning, INaab fc BurkerU (2003) argue that large ellip- 
ticals are formed by the dissipationless merger of two 
comparable-mass disk galaxies, i.e., major mergers, while 
low-luminosity ellipticals are produced by the dissipa- 
tionless merger of unequal mass disks, i.e., minor merg- 
ers. To demonstrate this possibility, INaab fc Burkerti 
( 2003) used numerical simulations of dissipationless disk 
galaxy mergers to show that remnants of major mergers 
rotate very little and are, in general, boxy, similar to lu- 
minous ellipticals. On the other hand, the simulated mi- 
nor mergers rotate significantly and have disky isophotes, 
similar to low-luminosity ellipticals. However, these sim- 
ulations did not include gas physics, and hence would not 
induce starbursts, quasar activi ty, nor satisfy the scalin g 
relations of elliptical galaxies ([Robertson et al. 1 12006a L 
Moreover, if most ellipticals are relatively old, i.e. have 
mainly old stellar populations, and were formed long ago 
by mergers, it is likely that the progenitor galaxies would 
have been gas-rich and had a higher gas fraction than lo- 
cal large, star-forming galaxies. 

In this paper we use a large suite of numerical simu- 
lations to explore the kinematic properties of remnants 
produced by the merger of comparable-mass disk galax- 
ies. We specifically address the differences between rem- 
nants formed via dissipationless mergers versus those 
produced in gas-rich mergers that include the cooling 
of gas, star formation and feedback. We find that gas- 
rich mergers can successfully reproduce many of the kine- 
matic properties of observed elliptical galaxies, while dis- 
sipationless remnants provide a poor match the data. 

The organization of the rest of this paper is as follows. 
In § we describe the numerical simulations, including 
the disk galaxy models and the galaxy collisions i|2.1|) . 
followed by the techniques employed to analyze individ- 
ual merger remnants l|2.2|l . In §[3 we present the results 
of our analysis for the entire series of merger remnants. 
To begin, we report the aggregate of measured properties 



for our entire series of simulated merger remnants i|3.1[l . 
Following this, wc show the remnant rotational support 
(|3.2|l . including where this rotation originates (|3.3I) . and 
what individual remnants are like l|3.4|l . To better char- 
acterize the remnants we also analyze their anisotropy 
(jEU, shape (E3J), minor-axis kinematics (|3.7(l . and how 
these quantities depend on our various input assump- 
tions H3.9|) . Finally in §0]we discuss the implications of 
our results for the formation of elliptical galaxies and we 
conclude in §0 

2. METHODS 
2.1. Merger Simulations 

The simulations presented here are part of a large, 
ongoing effort to study galaxy mergers and how this 
process impacts the formation and evolution of galax- 
ies. Hence, the methods used here are identi- 
cal to, and also described in, several related works 
(iCox et all 120061: iDi Matteo et all 120051 iHopkins et al] 
|2005al I2006al l2f)05bld l2006blcL l2005dt [Lidz et alJl200fit 
[ Robertson et alJ l2006al 120051 l2006bl: iSpringel et alJ 
2005aEjr Readers that desire a more detailed descrip- 
tion of the simulation code and the cons truction of the 
progeni tor galaxy models are referred to ISpringel et alJ 
(2005b), where the majority of the methods were first 
introduced. Below, we provide a brief overview of our 
methodology and focus on the assumptions most rele- 
vant for the results reported here. 

The simulations described here were performed us- 
ing the N-body/SPH (s moothed particle hydrodynamics) 
code GADGET2 l[Springell2005l) . w hich is based on a fully 
conse rvative formulation of SPH ( S pringel &: Hernauistl 
2002). Along with the standard features of this pub- 
lically available code, the version we employ includes 
radiative cooling of gas, star formation th at is de- 
signed to match the observed Schmidt-law l) Schmidt! 
1959; Kcn niguja , t he m ultiphase feedback model 
of Springel & Hcrnquist (2003) softened (qEOS=0-25) so 
that the mass- weighted interstellar medium temperature 
is ~ 10 4 ' 5 K, and a centrally located sink particle that 
can accrete gas and release isotr opic thermal energy th at 
represents a massive black hole l|Springel et al.ll20 05b). 

The fiducial galaxy model consists of a dark matter 
halo and an embedded rotationally-supported exponen- 
tial disk. For simplicity we do not include a spheroidal 
bulge in our fiducial galaxy model. T he dark matter halo 
is initialized with a lHernauisti l[1990t) profile that has an 
effective concentration of 9, a spin parameter A = 0.033, 
and a circular velocity V200 = 160 km s -1 . The expo- 
nential disk size is fixed by requiring that the disk mass 
fraction (4.1% of the total mass) is equal to the disk 
spin fraction. This results in a radial disk scale length 
of 3.9 kpc. In addition, a specified fraction / of the 
disk mass is in a gaseous component. The galaxy models 
are realized by 120,000 dark matter particles, and 80,000 
disk particles. A fraction / x 80, 000 of the disk particles 
represent the gaseous disk, and the remainder represent 
the collisionless stellar disk. 

In the present study, we restrict our analysis to simu- 
lations of major mergers between identical disk galaxies. 
We also consider only two types of progenitor disks. One, 
where the exponential disk is purely stellar (/ = 0), and 
a second where the disk is 40% gas (/ = 0.4). In all 
simulations, we adopt a gravitational softening length 
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merger simulations. Col. (1) is our unique orientation 
identification. Cols. (2) and (3) are the initial orientation 
of disk 1 and cols. (4) and (5) are the orientation of disk 2. 
A brief description of several unique orientations is listed 
in column (6) . 



of 140 pc. We restrict our analysis to quantities mea- 
sured on the scales of the effective radius, i.e., several 
kpc, where we are confident that we have sufficient res- 
olution. Unfortunately, with this resolution we cannot 
reliably determine the inne r (< 140 pc) surface density 
profiles l)Lauer et al.l 120051 and references therein) nor 
follow the merger of the bin ary black holes (see, e.g., 
Milo savlievic fc Merrittll200lD . 

Once built, pairs of identical galaxies are placed on 
parabolic orbits with the spin axis of each disk specified 
by the angles 9, and <fi in standard spherical coordinates. 
Table 2] lists the adopted orientation of each progeni- 
tor for our fifteen simulated mergers. The orientation 
list contains seven idealized mergers (labeled b — h) that 
represent orientations often seen in the literature. For 
example, orientation h has the galactic spin and orbital 
angular momentum aligned in what is typically called a 
prograde-prograde merger. Also included are slight vari- 
ations of this case, e.g., prograde-retrograde, retrograde- 
retrograde, and pola r. The remai nder of the mergers 
(labeled i — 6) follow iBarnesI l)1992|) in that they are se- 
lected to be unbiased initial disk orientations according 
to the coordinates of two oppositely directed tetrahe- 
dra. These orbit s are id entical to those explored by e.g. 
INaab fc Burkertj ((2003) and therefore allow for a direct 
comparison to their results. 

The primary results of our work are based upon two se- 
ries of simulated major mergers. The first series consists 
of the fifteen orientations listed in Table ^ each simu- 
lated with progenitor disks composed of purely collision- 
less stellar mass (/ = 0). The second series is identical to 
this first series but the progenitor disks all contain 40% 
gas (/ = 0.4). In what follows we will consistently refer 
to the former series as "dissipationless" mergers, while 
the latter are "dissipational" , or "40% gas" mergers. 

In § 13.91 we consider additional simulations to address 
how robust our results are to the chosen initial condi- 
tions. 



2.2. Remnant Analysis 

The techniques we employ to analyze the merger rem- 
nants are designed to mirror those typically adopted by 
observers. For each remnant we project the stars onto 
a plane as if observed from a particular direction. An 
example of the projected stellar mass is shown in the 
upper-left plot of Figure^ The stellar mass is composed 
of two components. One, the dissipationless disk parti- 
cles that are present in our initial conditions, and two, 
stars that are formed during the simulation from dense 
gas, when included. 

Once the stellar mass is projected onto a plane, we 
determine the iso-density contour that encloses half of 
the stellar mass and fit an ellipse to it. We label the semi- 
major axis of this ellipse a, and the semi-minor axis b. 
The ellipticity is then defined as e = 1— b/a. We note that 
a is slightly larger than the half-mass radius R e , which is 
typically calculated using a circular aperture. Both the 
circular-aperture half-mass radius, and the fitted half- 
mass ellipse are overlaid on the projected surface density 
in Figure ^ in addition, the iso-density contours are 
shown. The remnant shown in Figure ^ is a dissipational 
merger from the / orbit. It is a fast rotator viewed from 
an angle perpendicular to the orbital plane. 

We note that observations often determine R e in a 
manner which differs from our approach. In that case, 
the surface brightness profile is fit to an analytic function, 
e.g., a i? 1 / 4 ijde Vaucouleursjl94c| or Sersic ijSersid 1968(1 
profile, from which the effective radius is extracted. Un- 
fortunately, owing to the limited resolution of telescopes 
and finite backgrounds, this technique is subject to un- 
certainties that depend on the radial range over which the 
fit is performed and the assumed pro f ile ijTruiillo et all 
12001^) . I n fact. IBoylan Kolchin et all <|2005() show that 
the R e of lHernauistl 1)19 90) profile can be underestimated 
by nearly 30% when the maximum radius used for the 
fitting produces is varied between ~ 2 R e to > 10 R e . 
Because of these complications, and because we know 
the exact location of all stellar material we can extract 
the half-mass radius directly from the particle informa- 
tion. This eliminates any ambiguity between the pro- 
cedure used to fit the profile and the resultant R e . In 
what follows we will typically use a, rather than i? e , as 
a measure of the "size" of our merger remnants. 

To quantify the kinematics of each remnant, we place 
a slit along the major axis and measure the velocity and 
dispersion along the slit. Again, this is demonstrated in 
Figure ^ where we show the projected two dimensional 
rotation and dispersion velocity fields, the slit, and the 
resultant velocity profiles. The slit has length and width 
of ~ 3a, and 0.25a, respectively. The slit is divided into 
26 bins, lengthwise, and the line of sight velocity field is 
measured for each bin. For simplicity, the velocity profile 
in each bin is assumed to be Gaussian, and thus we ex- 
tract the mean rotation velocity and velocity dispersion 
from each bin. 

The rotation velocity U ma j along the major axis is de- 
termined as the average absolute value of the maximum 
and minimum velocity along the slit. The velocity dis- 
persion a of each remnant is the average of all dispersions 
R < 0.5a along the slit. In practice, the dispersion pro- 
file is fairly fiat out to ~ R e and thus the choice of 0.5 
does not affect our results. The large aperture used to 
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Fig. 1. — Projected mass density (upper left) and velocity (upper right) and dispersion (lower right) fields for one of our merger remnants 
(dissipational merger, orientation /, viewed perpendicular to the orbital plane) with overlaid isodensity contours. The half-mass radius, 
computed within a circular aperture, is shown by a (blue) circle. The half-mass isodensity contour has been fit with a (red) ellipse. The 
semi-major axis length a and ellipticity e are shown in the upper-left of the figure. A slit has been placed along the major axis and the 
velocity along the slit is show in the lower left panel. This remnant is a fast rotator. 



measure a is also selected becuase it is much greater than 
the numerical resolution. 

We note that the two dimensional velocity helds con- 
tain a wealth of information that is largely lost by the 
use of a simple slit. However, the majority of velocity 
data taken to date has been obtained using slits and we 
therefore follow this procedure in our work here. We 
do note, though, that a fruitful avenue for future work 
is the comparison of simulated two dimensional velocity 
fields such as those present ed in FigurelHto observ ational 
samples such as SAURON l)Car)oellari et al.l200filK which 
use integral field spectroscopy and thus image the two- 
dimensional velocity field. 

Because our models do not account for cosmological 
evolution, the simulation time is not directly connected 
to redshift. The initial conditions are designed to ap- 
proximate our own Milky Way galaxy except, when in- 
cluded, a higher disk gas fraction. Thus, the mergers 
presented here are probably more representative of mod- 
erate redshifts where disks were likely gas-rich. In fact, if 
the Milky Way has been forming stars at a steady rate of 
lM0yr~ 1 , without any other mass flux, it would have had 
a gas fraction above 40% about 10 Gyr ago (z ~ 1.5). In 



any case, we will mainly compare our remnants to local 
elliptical galaxies, most of which are old stellar systems 
and dynamically relaxed. Hence, we would like to ensure 
that the merger remnants are relaxed systems as well. In 
order to investigate the dynamical stability of our rem- 
nants, Figure [3 shows the size, both a and R e , ellipticity 
e, major axis rotation speed Vmaj, central velocity dis- 
persion tr, and the ratio Vmaj/c as a function of time 
after the merger is complete for the same remnant that 
was shown in Figure ^ For the most part, the remnant 
relaxes quickly (< 300 Myr) after the merger. Prior to 
this, the central dispersion a, the rotation speed, and size 
appear to be the most strongly time-dependent quanti- 
ties. In particular, the dispersion is about ~15% larger 
immediately after the merger compared to the relaxed 
remnant. Of all the quantities, the rotation speed fluc- 
tuates the most, and it varies within 10% of its average 
value, throughout the simulation. These 10% fluctua- 
tions are carried over to the ratio Vmaj/c, and thus we 
consider this a reasonable quantification of the error as- 
sociated with when we observe each remnant, which we 
take to be ~ 2.5 Gyr after the merger is complete. 
As a final comment, we note that all results presented 
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remnant displayed in Figure 111 



in this work use the stellar mass as opposed to the stellar 
light. This choice is motivated in part by the uncertain 
star formation history of all particles designated as stars 
at the beginning of the simulation. As we show in § 13.31 
the stars formed during the simulation tend to rotate 
more than those present at the beginning of the simula- 
tion, and thus our choice to weight the velocity by mass 
as opposed to luminosity may underestimate the veloci- 
ties as compared to recent merger remnants while being 
a fair representation of most old ellipticals. 

3. RESULTS 

3.1. Effective Radius, Ellipticity, Dispersion, and 
Rotation 

Figure shows histograms of the semi-major axis a, 
half-mass isophote ellipticity e, central velocity disper- 
sion cr, and maximum rotation speed along the major axis 
Vmaj. Each of the 30 remnants (15 dissipational, and 15 
dissipationless) is projected along 195 lines of sight. The 
angle s are selected using the HEALPIX (G6rs kT"et alJ 
2005) software with nside = 4 and uniformly sample 
a unit sphere. This procedure results in a total of 2250 
"data" points which are displayed in FigureElfor both the 
dissipationless and the dissipational merger remnants, 
represented by the (blue) open and (red) cross-hatched 
histograms, respectively. 

It is apparent that dissipational merger remnants are 
quite different from their dissipationless counterparts. 



The effects of gaseous dissipation and star formation re- 
sult in remnants that are more compact, rounder, have 
higher velocity dispersion, and have a much more uni- 
form distribution of rotation speeds. 

The size of the remnants is a strong function of dis- 
sipation. In fact, every dissipational remnant is smaller 
(by nearly ~ 50%) than its corresponding dissipation- 
less version. This fact is unsurprising given the signif- 
icant gas fraction of our progenitor disks and the effi- 
cienc y with which mergers drive gas in to the galaxy cen- 
ters l|Barnes fc HernauistJ 119911 11996ft, where it f uels a 
burst of star formation IjMihos &: Hernauistj H996K For 
the mergers we present here, this process is so efficient 
that the average half-mass radius R e of the dissipational 
remnants, 4.0 kpc, is smaller than the original progeni- 
tor disk 4.4 kpc. The average R e of the dissipationless 
remnants is 6.9 kpc. The average semi-major axis a of 
the dissipational remnants is 4.1 kpc and 7.9 kpc for the 
dissipationless remnants. 

To check the absolute size of our merger remnants, we 
note that the size-mass relationship found in the Sloan 
Digital Sky Survey (SDSS) for early-type galaxies sug- 
gests a typical size of 3.4 kpc for galaxies of equivalent 
stellar mass (~ 8 x 1O 1O M0) to our simulated remnants 
(jShen et al.ll2003|) . Hence, it appears that both our dis- 
sipational and dissipati onless remnants are larg er than 
SDSS galaxies (see als o iRob ertson et al.ll20 Q6alK H ow- 
ever, as pointed out by Bovlan-Kolchin et al. ( 2005, Ap- 
pendix A), sizes can be systematically underestimated 
depending on the range of radii considered when fit- 
ting an assumed profile. We also note our initial disk 
galaxies are relatively large, a reflection of an above av- 
erage spin parameter A assumed for the dark matter 
halo. In reality there exists a distribution of spins and 
thus disk sizes. Furthermore, IRobertson et all ((2006a) 
have shown that the sizes of the remnants depend some- 
what on the orbit o f the merger, in agreement with 
iBovlan-Kolchin et alJ l)2005|) . and matching the obser- 
vations would thus require a distribution of orbits con- 
sistent with those expected cosmologically. Therefore, 
we consider this (dis)agreement as marginal at best and 
here we simply note that the dissipational remnants are 
in closer agreement with the SDSS size-mass relation for 
early- type galaxies than dissipationless remnants. 

The top-right plot in Figure [31 shows the distribution 
of ellipticities for our merger remnants. Ov erplotted are 
150 ell ipticities as measured and compiled bv lFranx et alJ 
i!991t). and anoth er 420 elliptical galaxies observed by 
Fabe r et al.l Jl989). The dissipational merger remnants 
are nearly indistinguishable from the two samples of ob- 
served ellipticals. Both the observed and dissipational 
ellipticity distributions peak near 0.2 and then drop pre- 
cipitously above 0.4. There exist almost no highly el- 
liptical (e > 0.6) galaxies. These features are in sharp 
contrast to the dissipationless remnants whose distribu- 
tion of ellipticities peaks near 0.4 and exten ds bey ond 
0.7, similar to that found bv lNaab fc BurkertJ <j2003|h 

The central velocity dispersion is a common measure 
of galaxy mass and appears to be correlated with many 
fundamental properties of elliptical galaxies. As the his- 
togram of remnant velocity dispersions in Figure|3shows, 
the central velocity dispersion is highly dependent upon 
viewing angle, merger orbit, and the dissipative effects 
of gas. The distribution of dissipationless dispersions is 
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Fig. 3. — Histograms of remnant properties. Shown in (red) cross-hatch are the dissipational remnants, while the open (blue) histogram is 
for the dissipationless remnants. Histograms show, clockwise from upper-left, the semi-major axis a, the half-mass isophote ellipticity e, the 
central velocity dispersion a, and the maximum velocity along the major axis V ma j. Each histogram represents the composite properties 
of all fifteen remnants, viewed from 150 projections. The semi-majo r axis histog r am in clud es the half-mass rad ius of the progenitor disk. 
The ellipticity histogram includes data on observed ellipticities from Fa ber et alJ I198SA and|Franx ct al. ( 1991J). 



well fitted by a Gaussian with mean 150 km s _1 and vari- 
ance 11 km s . The distribution of dissipational merger 
remnants is better described as bi-model, with one com- 
ponent a Gaussian centered at 178 km s _1 and variance 
10 km s _1 , and a second component that extends to 
lower velocity dispersion. This second, lower dispersion, 
distribution is comprised of five mergers (d,f,i,l,m) that 
are all significant rotators. The mean and variance of 
the entire dissipational distribution is 169 km s and 
15 km s _1 , respectively. 

Because it probes the gravitational potential, the cen- 
tral velocity dispersion is often used in combination with 
the effective radius as a proxy for galaxy mass (i.e., 
M oc a 2 R e ). While it does not affect the results of this 
paper, we note that a 2 R e for our remnants can vary by 
a factor of 2, even though all the remnants galaxies have 



exactly the same total mass. 

The increased velocity dispersion of the dissipational 
remnants is accompanied by a steeper surface density 
profile, as shown in Figure 0| Within a radius of ~ 
1.8 kpc the dissipational remnants have a much higher 
surface density owing to gas dissipation and star forma- 
tion. This feature is present regardless of the merger or- 
bit, even though the different orbits produce a variety of 
central (< 1 kpc) profiles, ft is intriquing that the surface 
densities appear to be well fit by a i? 1 / 4 -profile and do 
not demonstrate the distinct cusp of newly formed stars 
seen in previous simulations |Mihos fe HernauisH Il994at 
ISnrmgell200fl ICox et al] 12005ft . Future work is under- 
way to address whether these surface density profiles are 
consisten t with observed ellipticals and merger remna nts 
(see, e.g.. lRothberg fc Joseol]l2004HLauer et alJl2005D as 
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Fig. 4. — The surface density profile for all thirty (the fifteen 
orientations from Table for both the dissipational and dissipa- 
tionless mergers) remnants plotted against il 1 / 4 as viewed from a 
single projection normal to the orbital plane. The dark solid lines 
that have uniformly higher central surface density are the dissipa- 
tional merger remnants. The light dashed lines with lower central 
surface density of the dissipationless merger remnants. The dashed 
vertical line denotes the softening length, and hence the resolution 
of the simulation. 



well as whether the profiles depend on progenitor disk 
gas fraction or mass. 

Spheroidal (and possibly all) galaxies are thought to 
host black h oles in their centers in accor d with the Mrh — 
a re lation (jFerrarese h Merrittl 1200(1 iGebhardt et all 
2000). For the remnants displayed here, with central 
velocity dispersions ^150 km s , each remnant should 
contain a fairly massive (~ several xl0 7 M Q ) black hole. 
During the dissipational mergers, black holes grow to 
the appropriate mass and futher accretion is terminated 
by a feedback induced galactic wind. Our merger are 
thus consistent with the idea that the Mbh — cr re- 
lation is p roduced by the self-regulated growth of the 
black hole iDi Matteo et alJl2005t iSpringej et alJl2005ht 
iRobertson et alJ l2006tA iKa zantzidi s et al.ll2005|) . Be- 
cause of the lack of gas, the black holes in the dissipa- 
tionless mergers remain the same size during the entire 
simulation. In this case, the initial black hole masses 
require some fine-tuning in order for the remnants to re- 
side on the observed Mbh — cr relation. We consider this 
further evidence for the necessity of gas-rich progenitor 
disk galaxies. 

The final histogram in Figure |3] shows the maximum 
rotation velocity along the major axis of each remnant. 
As with the previous histograms, the dissipational and 
dissipationless remnants are quite different. While the 
dissipationless remnants rotate very slowly, with all mea- 
surements below 50 km s _1 , the dissipational remnants 
span a wider range of rotation speeds, with some projec- 
tions reaching speeds of ~ 140 km s _1 . We will discuss 
the rotational properties of our merger remnants in more 
detail in the following section. 

3.2. Rotational Support 



As stated in the introduction, elliptical galaxies can 
be split into two groups. Small ellipticals tend to ro- 
tate along their major axes, while large elliptical galaxies 
show little or no rotation. More specifically, the rotation 
of small ellipticals is consistent with that of an oblate, 
isotropic rotator. Thus, their elliptical shape can be un- 
derstood solely i n terms of flattening induced by their ob - 
served rotation <)Binnevll978HBinnev k. Tremainell987|) . 
Large ellipticals, on the other hand, have little rotation. 
Their elliptical structure must be supported by some- 
thing other than rotation, and is most likely velocity 
anisotropy. One of the clearest ways to visualize this 
trend is to plot the rotation velocity Knaj divided by the 
central velocity dispersion a, against the ellipticity e, in 
what is commonly referred to as the anisotropy diagram. 

Figure El shows our merger remnants in an anisotropy 
diagram. For reference, overplotted on this dia- 
gram are data from several observational samples 
lIDavies et alJll983t IBenderlll98l iBender fe Nietolll99tt 
Ide Zeeuw et al.H2002|) . IBinnevI l|1978|) has shown that for 
oblate isotropic rotators V/a is purely a function of ellip- 
ticity e. If this is the case, then at high ellipticities one is 
viewing a rotating, flattened isotropic ellipsoid from edge 
on. The solid line in Figure is this relation for oblate 
isotropic rotators 

(V/a) = [e/(l-e)] 1 / 2 . (1) 

ijBinnevI Il97cf) . In terms of the anisotropy diagram, 
galaxies that are near the line are presumably oblate 
isotropic rotators. Their shapes are flattened by rota- 
tion and they are being viewed from a variety of angles 
with respect to their rotation axes. Galaxies that are 
below this line, i.e. those that have very little rotation, 
and some amount of ellipticity, arc likely to have strong 
velocity anisotropics. 

As in the previous section, the dissipationless and dis- 
sipational remnants occupy very different regions in the 
anisotropy diagram. The dissipationless remnants, pre- 
sented in the right-hand plot of Figure03 occupy a nearly 
identical region of the anisotropy diagram as those pre- 
sented in lNaab fc BurkerU l)2003|) . Consistent with their 
findings, most r emnants are very elliptical and hardly ro- 
tating (see alsolNegroponte fc Whitelll983l !Barneslll983 
iHernauistl I1992T . We note that our additional orbits 
(b — h) produce non-rotating remnants exactly like orbits 
i — p ( i.e., the exact orbits simulated bv lNaab fc BurkertJ 
(2003)). There are a large number of remnants that 
have an elliptici t y grea ter than 0.6, as was noted by 
iNaab fc Burkertl pOO^ . and was clearly displayed in 
the histogram of Figure |3| This distribution of ellip- 
ticities is not observed in real ellipticals, especially for 
galaxies with little or no rotation. Ellipticity aside, the 
small rotation of the di ssipationless merger remnants led 
INaab fc Burkertl (|2003) to suggest dissipationless merg- 
ers as a possible formation mechanism for bright ellipti- 
cals. 

The dissipational merger remnants are shown in the 
anisotropy diagram on the left panel of Figure These 
remnants span a much larger range of rotational prop- 
erties than the dissipational remnants, and include rem- 
nants that reside on the oblate isotropic rotator line as 
well as remnants that fall below this line. The most com- 
mon position for dissipational remnants to reside in the 
anisotropy diagram is V/a ~ 0.35 and e ~ 0.18. At 
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Fig. 5. — The V ma j /or versus ellipticity diagram for dissipational (40% gas) and dissipationless merger remnants. V ma j is the maximum 
rotation speed measured in a slit along the major axis, a is the velocity dispersion a vera ged within half of an half-mass radius, and the 
ellipticity is measured at the h alf-mass isoph ote. Further details can be found in § 12.21 The soli d line in both plots i s that expected 
for an oblate isotropic rotator IBinneyl 119781). Overplotted arc data from observed ellipticals from [Davics et alj Jl98lI) ; IEenderl (1988); 
IBender fc Nietol tl990h : IdeZeeuw et alj 1200211 . 
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Fig. 6. — Histogram of (V ma j/o")* for both the dissipationless 
(open, blue) and dissipational (filled, red) merger remnants. Also 
shown is a histogram of observed ellipticals for the same points 
shown in F igure 151 and a histog ram of the sample of spheroids 
compiled bv [Bender et all ^992). 



this point the remnant is slowly rotating, is below the 
oblate isotropic rotator line, and its structure is presum- 
ably supported by anisotropy; all of these are features 
similar to bright ellipticals. 

In general, Figure shows that the dissipational rem- 
nants encompass a significant fraction of the observa- 
tional data points; much more so than the dissipation- 
less remnants. To investigate this further, for each 
one of our merger remnants we evaluate (Vmaj/er)*, the 
ratio of the measured rotation parameter V ma j/a to 
Equation QJ, the value for an isotropic oblate spheroid 



flattened by rotation (iBinnevI 119781: iDavies et"aTlll983t 
IKormendv & Benderl H996) . In a manner similar to Fig- 
ureEJwe show histograms of (V ma j/a)* for both the dissi- 
pationless and dissipational merger remnants in Figure^ 
In addition we show the distribution for the same obser- 
vations as shown in Figure |S] as well as the observational 
sample of IBender et al.l l|1992j) . We caution that these 
samples are not independent of one another as there are 
a number of galaxies which are counted in both. Con- 
firming the trend in Figure we see that the observa- 
tions and dissipational remnants have a predominantly 
flat distribution below (V ma j/a)* < 1, and a tail of val- 
ues above this. In contrast, the dissipational remnants 
are peaked below 0.2 and there are almost no remnants 
above 0.6. 

Even though the dissipational remnants are a 
much better match to the observed distribution of 
(Knaj/c)*, as shown in Figure El there is still a dis- 
crepancy between the observed and predicted num- 
ber of fast rotators; i.e., those with (V ma j/cr)* > 
1. This deficiency was also noted by several other 
studies under a v a riety of different circumstances 
llCretton et all [200lt IGonzalez-Garria k Balcellsl l2005t 
I Gonzalez-Garcia fc van Albadal 120051 INaab fc BurkertJ 
2003), even when minor mergers were included. Observa- 
tions find a large number of galaxies exhibiting fast rota- 
tion. In the observations plotted here, 20% of the galax- 
ies shown in Figure have ( V m ?,;JcrY > E and th is num- 
ber increases to 33% for the Bender et al.l ()1992|) sample 
which includes a larger number of faint ellipticals and 
bulges. The fraction of fast rotators could increase even 
more with the inclusion of SOs and o bservations whic h 
probe the velocity field at large radii (|R.ix et al.l tl999). 
For our simulated dissipational merger remnants, only 
11% of the projected images have (V ma j/cr)* > 1. For 
the moment, we do not find this discrepancy alarming 
as we have sampled a limited range of disk galaxy initial 
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Fig. 7. — The rotation curve for the dissipational merger rem- 
nant /, viewed perpendicular to the orbital plane. This is one of 
the fastest rotating remnants and is identical to that used in Fig- 
ures [I]and|2] The rotation curve is decomposed by particle type 
according to the key inset within the bottom-left of the plot. The 
dissipational stellar particles are segregated based upon their rela- 
tive age; "old" stars begin the simulation as a collisionless particle, 
"intermediate" stars were created during the isolated and first pas- 
sage phases, while "young" stars are created from dense gas during 
the final merger. For reference the composite "all stars" shows the 
total stellar rotation curve, and the solid black line shows the dis- 
sipationless merger remnant viewed from the same angle with an 
identical slit. 



conditions. In particular, we suspect that varying the 
gas fraction may produce a much larger fraction of fast 
rotating remnants, continuing the trend seen in Figure^] 
We briefly address this in § 13.91 

3.3. The Origin of Stellar Rotation 

To understand what generates the significant rotation 
found in dissipational remnants, we decompose the stel- 
lar rotation curve presented in the lower-left of Fig- 
ure ^ into contributions from stars by their age. The 
resulting rotation curves are shown in Figure {7\ Star 
particles are designated as "old", "intermediate", or 
"young" on the basis of when they originally became 
a dissipationless stellar particle. Old stars are those 
that are present prior to the start of the simulations, 
and thus designated as stars when setting up the ini- 
tial conditions. Star formation during the course of the 
merger simulation predominantly occurs in bursts that 
accompany the first passag e of the galaxies, and the 
final merg er eve nt (see e.g.. |Miho^^ ^Hernauisti ri994c: 
iCoxet all 120051: iSnringel et all l2005bD . These burst 
events occur ~1 Gyr apart. We denote any stars created 
during the isolated evolution prior to the first burst, and 
those generated during the first burst as "intermediate" . 
All stars created during the final merger are "young" . In 
addition to the rotation curves for stars in the dissipa- 
tional run, we also plot the stars from the dissipationless 
simulation as a solid black line in Figured 

Figure [7| demonstrates that it is the stars born dur- 
ing the final merger, i.e., the "young" stars, that are the 
fastest rotators. Because this young stellar component is 
centrally concentrated (its half-mass radius is less than 
1 kpc) and is the dominant stellar component inside of 



~1 kpc, the central rotation closely tracks that of the 
young stars. Outside of 1 kpc the "intermediate" and 
"old" populations dominate the stellar mass and thus 
the rotation curve closely tracks these components. Both 
"intermediate" and "old" components have similar rota- 
tion suggesting that all stars that exist prior to the final 
merger are equivalent in terms of their rotation proper- 
ties in the merger remnant. 

Figure also plots the rotation curve for the dissipa- 
tionless version of this merger, shown with a solid black 
line. While the dissipationless remnant also exhibits ro- 
tation, the rotation curve is much shallower than the dis- 
sipational remnant. Obviously, gaseous dissipation and 
star formation have affected the rotation of the old stellar 
population in the dissipational remnant. The most plau- 
sible explanation for this is simply the conservation of 
angular momentum. The central concentration of young 
stars steepens the potential well and draws the old stars 
into the center. The half-mass radius of the "old" stars 
is ^20% less than in the dissipational remnant and thus 
the velocity of these stars must be roughly ^20% higher 
to conserve angular momentum. The fact that the veloc- 
ity increase is larger than this may represent a transfer 
of angula r momentum between the s tars and gas, as was 
shown bv lMihos fc HernauistJ 1^)96). 

Even though the young stars rotate faster than the 
old stars, most of the angular momentum is carried by 
the old stars. Specifically, the old stars typically have 
three times the specific angular momentum of the young 
stars. However, most of the angular momentum carried 
by the old stars is at radii larger than the effective ra- 
dius, and thus is not contained in our slit which only 
extends slightly beyond one effective radius. Figure 
shows that beyond ~4 kpc, the dissipational and dissipa- 
tionless remnants have similar rotation, and also similar 
angular momentum distributions for the old stellar com- 
ponent. Thus, these results suggest that the significant 
rotation seen in dissipational remnants is not indicative 
of more angular momentum, but dissipation and star for- 
mation effectively redistribute angular momentum from 
large radii to within the effective radius where rotation 
is typically measured. 

Although our results indicate a segregation between 
stellar age and rotation it is unclear if this trend could 
be verified by observations. As the stellar popula- 
tions evolve, the small age differences (1-2 Gyr) will 
quickly be washed out and thus difficult to disen- 
tangle. However, the general trend for young stars 
to be centrally concentrated and kinematically differ- 
ent than the more extended older stars is a generic 
feature of these simulations. In a few cases the 
merger results in counter-rotating or kinematically de- 
coupled cores, as is seen in many observed e llipticals 
Ipavies fc Birkirishawlll988t iFranx et al.lll989j) . and as 
found earlier bv lHernauist fc Barnesf ifl99Fl using simu- 
lations with a smaller gas fraction than here and ignoring 
star formation. 

3.4. Individual Merger Remnants 



Although dissipational remnants display much more 
rotation, on average, than their dissipationless counter- 
parts, there is a large spread of rotation properties be- 
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TABLE 2 

Projection Averaged Remnant Properties 



Orientation I Dissipationless 40% Gas 



ID 


c 


Re 


a 




cr 


(V maj /<x)* 


e 


Re 


a 


Vnaj 


a 


0W<t)* 






(kpc) 


(kpc) 


(km s" 1 ) 


(km s -1 ) 






(kpc) 


(kpc) 


(km S" 1 ) 


(km s —1 ) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


h 


0.48 


6.49 


8.73 


14.3 


150.8 


0.09 


0.27 


3.56 


4.16 


46.0 


180.7 


0.25 


b 


0.46 


6.27 


8.21 


11.7 


150.0 


0.08 


0.27 


3.70 


4.21 


26.8 


182.2 


0.15 


c 


0.50 


6.37 


8.39 


15.8 


153.5 


0.10 


0.32 


3.70 


4.43 


116.6 


177.5 


0.65 


d 


0.42 


6.40 


8.04 


9.1 


151.0 


0.06 


0.24 


3.86 


4.34 


111.9 


147.7 


0.76 


e 


0.37 


6.23 


7.87 


9.7 


149.7 


0.06 


0.21 


3.49 


3.99 


11.3 


178.5 


0.06 


f 


0.33 


6.16 


7.17 


15.3 


147.0 


0.10 


0.26 


3.99 


4.54 


143.3 


140.5 


1.03 


g 


0.37 


6.19 


7.60 


13.0 


150.2 


0.09 


0.19 


3.84 


4.19 


65.5 


172.7 


0.37 


i 


0.34 


6.24 


7.71 


13.8 


150.4 


0.09 


0.19 


3.64 


4.00 


115.7 


155.4 


0.75 


j 


0.22 


6.09 


6.91 


15.5 


149.6 


0.10 


0.13 


3.36 


3.63 


66.0 


169.6 


0.39 


k 


0.49 


6.39 


8.83 


9.3 


151.4 


0.06 


0.22 


3.56 


4.01 


73.2 


175.2 


0.41 


1 


0.30 


6.24 


7.49 


27.2 


147.8 


0.18 


0.19 


4.12 


4.49 


91.5 


147.8 


0.62 


m 


0.39 


6.21 


7.99 


14.0 


151.4 


0.09 


0.17 


4.01 


4.30 


77.2 


146.2 


0.53 


n 


0.26 


6.29 


7.34 


12.3 


150.4 


0.08 


0.21 


3.50 


3.93 


131.4 


169.4 


0.77 


o 


0.35 


6.20 


7.46 


31.2 


146.6 


0.21 


0.20 


3.67 


4.11 


128.7 


166.4 


0.77 


P 


0.45 


6.29 


8.21 


8.0 


152.1 


0.05 


0.24 


3.59 


4.07 


66.0 


168.6 


0.39 



Note. — List of properties of simulated merger remnants averaged from 190 viewing angles selected to uniformly sample the unit sphere. 
Column (1) is the ellipticity of the half-mass isophotc. Column (2) is the half-mass radius as measured using circular apertures. Column (3) is 
semi-major axis of the half-mass isophotc. Column (4) is the maximum rotation speed along the major axis. Column (5) is the central velocity 
dispersion, averaged within 0.5a, and column (6) is the ratio of V ma j/(7 to that expected from an oblate isotropic rotator. 



tween individual remnants ijBendo fc Bar nes 2000). Fig- 
ure [S] demonstrates this behavior by displaying the rota- 
tional properties of individual merger remnants, one per 
panel, viewed from 190 projections. It is these remnants 
that are co-added to produce Figure El 

Figure |S| clearly shows that projection effects cause 
each remnant to occupy a band in the anisotropy di- 
agram. Every remnant appears nearly circular from a 
small range of viewing angles and elliptical from most 
viewing angles. It is typically the case that remnants ro- 
tate very little when they appear circular and thus each 
band passes near the origin of the anisotropy diagram. 
At larger projected ellipticities there exists a wide range 
of rotational properties for the fifteen merger remnants. 
In some cases (d, /, i, and o) the remnant closely re- 
sembles an oblate isotropic rotator. However, in most 
cases the remnant swath is below the oblate isotropic ro- 
tator line. In all cases (except, perhaps, the j merger 
remnant) the maximum rotation is coincident with the 
largest projected ellipticity, suggesting that the majority 
of the dissipational merger remnants are (nearly) oblate 
and the rotation axis is along their short axis. 

It is important to emphasize the large range of kine- 
matic properties present in the dissipational merger rem- 
nants. This range of remnant properties owes entirely to 
the variations in initial disk orientations. Specifically, it 
appears that the co-planar orbits (h, b, and c) produce 
very flattened anisotropic remnants, while polar orbits, 
i.e., where one disk orbits in a 90° orientation to another 
disk (d and /), result in fast rotating remnants that are 
oblate and isotropic. Finally, tilted orbits, such as e, j, 
and k produce fairly round remnants that are not rotat- 
ing. 

The dissipationless merger remnants shown in Figure|Sl 
show much less variability than the dissipational rem- 
nants. In all cases, the remnants span a wide range of 
ellipticities and are never rotating. Elliptical remnants 
that do not rotate are thought to owe their ellipsoidal 



shape to anisotropic velocity dispersions. We will inves- 
tigate this in the following section. 

A final point concerns the initial disk orientations and 
whether or not they faithfully represent the full range of 
possibilities. Orientations i — p were originally selected 
to be a uniform sampling of disk orientations fsee lBarnesl 
1992). However, as we note above, the co-planar orbits 
(h, b, and c) and polar orbits (d and /) produce some- 
what peculiar remnants that are not well represented by 
the uniform sampling introduced by Barnes. Hence, it 
seems that one must include both a uniform grid to rep- 
resents the majority of merger events as well as a small 
sample of hand-selected orientations to represent more 
pathological orbits in order to gain a full understanding 
of the various properties possible in the aftermath of a 
galaxy collision. 

3.5. Anisotropy 

In the previous section, we demonstrated that all of 
the dissipationless and several of the dissipational merger 
remnants are ellipsoidal and slowly rotating. In short, 
they reside below the oblate isotropic rotator line in the 
anisotropy diagram. In this case, it is assumed that the 
ellipsoidal shapes are supported by an anisotropic ve- 
locity dispersion. However, it has recently been shown 
that mergers between unequal mass gal axies can produce 
fast r otators that are very anisotropic ijBurkert &: Naahl 
|2005|) . and thus we wonder if the slowly rotating rem- 
nants are necessarily anisotropic. In this section we at- 
tempt to address this question by measuring the velocity 
anisotropy of each merger remnant. 

We note that any dissipationless system in equilibrium 
can be accurately described by the tensor virial theorem 

2T jk + U jk = -W jk , (2) 

where Wjk is the potential energy tensor, and Tj k and 
n jk are the ordered and random components, respec- 
tively, of the kinetic energy tensor ijBinnev fc Tremaind 
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Fig. 8. — Vmaj/cr versus ellipticity, or anisotropy, diagram for individual merger remnants. The top grid of fifteen panels shows remnants 
where the progenitor disk contained 40% gas and the bottom grid is the corresponding mergers when the progenitor disks are dissipationless. 
Each plot shows one merger remnant for one of the merger orientations listed in Tabled viewed from 190 different viewing angles chosen 
to uniformly sample the unit sphere. The solid line i s tha t expected for an oblate isotropic rotator and the number in the middle-right of 
each panel is the anisotropy parameter 8 defined in § 13.51 For reference, <5 = is isotropic and 8 > is a sign of anisotropy. 



119871 eq. 4-78). Taking the trace of Equation J2J) re- 
sults in the scalar virial theorem, which simply relates 
the total kinetic energy to the total potential energy. 

To calculate the kinetic energy tensor we assume each 
N-body stellar particle has a gaussian velocity distribu- 



tion. The mean and distribution are computed by kernel 
weighting among the 96 nearest neighbors. A Cartesian 
coordinate system is selected such that the z-axis is par- 
allel to the original orbital angular momentum and the 
x and y axes are arbitrary. The mean velocity then con- 
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FlG. 9. — Histogram of remnant velocity anisotropy parameters 
5, where S is calculated via Equation J1J. With this definition for 
anisotropy, <5 = is isotropic and 8 > is anisotropic. 



tributes to Tjk and the dispersion to Ujk . For our merger 
remnants the kinetic energy is predominantly in the form 
ol random motions, with 



Trace(n jfc ) 
Trace(T jk ) 



> 3. 



(3) 



In comparison, this ratio for the progenitor disks is ~ 
0.17. 

The anisotropy of each remnant, then, is the degree to 
which the elements of ILjk, the velocity dispersion tensor, 
are unequal. Owing to our initial orbital configuration, 
in which the orbital angular momentum is aligned with 
the z-axis, we typically find H xx ,H yy > H zz . We follow 
iBurkert &: Naabl l)2005j) and define the anisotropy of the 
system as 

2TT 

8 = 1- = Z L . (4) 



n x 



With this definition, an isotropic system has (5 = 0, and 
8 > indicates some degree of anisotropy. For two of the 
dissipational remnants (g and n), 8 is slightly negative. 
However, these remnants are nearly isotropic and in this 
case we quote the absolute value of 8. 

A histogram of remnant anisotropics is presented in 
Figure El and the anisotropy of individual merger rem- 
nants is printed in each panel of Figure|Sl along with the 
swath of Vmaj/cr owing to various projections. From Fig- 
ure we see a systematic decrease in 6 for dissipational 
remnants compared to dissipationless remnants. In other 
words, every dissipational remnant is much closer to 
isotropic than its dissipationless counterpart. Roughly 
two-thirds (9/15) of the dissipational remnants are con- 
sistent with being isotropic (S < 0.05) compared to only 
one of the dissipationless remnants (J). Although there 
is a fair amount of scatter, there seems to be a uniform 
change in the anisotropy for dissipational mergers equiv- 
alent to AS ~ —0.16 compared to dissipationless rem- 
nants. 

As with the rotation, Figure[S]also shows that the ini- 
tial disk orientations have a large effect on the remnant 



anisotropy. This fact is perhaps not that surprising given 
the strong link between rotation and anisotropy. In par- 
ticular, the three initial orientations in which the initial 
disks are aligned with the orbital plane, i.e., co-planar 
mergers, (h, b, and c) result in remnants that are much 
more anisotropic than any of the other orientations. This 
trend appears to hold regardless of the presence of gas 
or not. 

Figure |H| also allows us to assess the relation between 
our (rough) measure of anisotropy 8 and rotation, as 
measured by the anisotropy diagram. As outlined in 
§ 13.21 these quantities should be anti-correlated, so that 
fast rotators are likely to be isotropic, while flattened 
non-rotators should be anisotropic. However, this rela- 
tion has not been tested in numerical experiments such 
as ours. 

At first glance, we note that the three dissipational 
remnants (d, /, and i) whose projected rotation is con- 
sistent with the line for oblate isotropic rotators (the 
solid line in the anisotropy diagram) are also isotropic 
(8 < 0.05). Next, many of the remnants that show 
very little projected rotation (dissipational remnants h, 
b, c, k, p and every dissipationless remnant except j) 
are all anisotropic. Therefore, for at least 22 of the 30 
total merger remnants, the expected relation between 
anisotropy and rotation appears to hold. 

However, there remain eight remnants for which their 
rotation properties are difficult to reconcile with their 
anisotropy. Of these eight, only one dissipationless rem- 
nant seems out of place; j. This merger remnant is not 
rotating and is moderately flattened (average cllipticity 
of 0.22) yet is apparently isotropic. One of the dissipa- 
tional remnants, e, is similar. The six remaining outlyers 
are all dissipational remnants (g, j, I, m, n, and o) and 
are all quite similar. These remnants are round (average 
ellipticity w 0.18), and have a large range of rotational 
properties for any one projected ellipticity. As such these 
remnants appear to occupy more extended regions within 
the anisotropy diagram, rather than distinct swathes. 
Future work will be necessary to determine whether the 
measure of anisotropy provided by equation Q is suffi- 
cient to describe these systems as well as to determine 
what generates their anomalous properties. 

3.6. Shapes 

The results presented in § 13.11 specifically the distri- 
bution of projected ellipticities in Figure [21 show that, 
in projection, dissipational remnants tend to be rounder 
than dissipationless remnants. To further quantify this 
trend we determine the three-dimensional shape of each 
rem nant as measured by the axis ratios of the inertia ten - 
sor ljBarneslll992l IHernauistllT99l Il99l ISpringell [2000). 
Because our analysis has focused upon the stellar com- 
ponent, primarily within R e , we use only the most bound 
half of all stellar particles to calculate the inertia tensor. 
We note that there are many other methods by which 
the shapes can be quantified. Most of these methods 
calculate the eigenvectors from some form of the iner- 
tia tensor. In our case, we calculated the inertia tensor 
based upon the most bound stellar particles but alter- 
natives are spherical or ellipsoidal windows which may 
be interatively adjusted until t he measured shape and 
selected aperture agree (see e.g- IAlkood et al.ll2006l for 
survey of different methods and recent results from the 
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TABLE 3 

Remnant Shapes 





Dissipationless 




40% Gas 


ID 


6 


c 


T 


b 


c 


T 


h 


0.73 


0.42 


0.56 


0.94 


0.43 


0.14 




0.86 


0.39 


0.30 


0.94 


0.42 


0.14 




0.75 


0.37 


0.50 


0.98 


0.40 


0.04 


d 


0.77 


0.55 


0.58 


0.90 


0.67 


0.34 


e 


0.76 


0.62 


0.68 


0.85 


0.71 


0.57 


f 


0.78 


0.67 


0.67 


0.94 


0.40 


0.41 


g 


0.79 


0.62 


0.61 


0.86 


0.74 


0.57 


i 


0.82 


0.60 


0.52 


0.93 


0.69 


0.24 


j 


0.84 


0.76 


0.73 


0.96 


0.82 


0.24 


k 


0.69 


0.49 


0.69 


0.93 


0.62 


0.21 


1 


0.91 


0.61 


0.27 


0.93 


0.71 


0.28 


m 


0.72 


0.60 


0.75 


0.88 


0.71 


0.44 


n 


0.81 


0.73 


0.74 


0.97 


0.81 


0.15 


o 


0.96 


0.62 


0.12 


0.96 


0.79 


0.19 


P 


0.76 


0.52 


0.58 


0.93 


0.63 


0.23 



Note. — List of remnant axial ratios b and c, 
and triaxiality parameter T, as defined in 



to assess. From the distribution of observed axial ra- 
tios, several works have shown that elliptical galax- 
ies can not be enti rely c o mposed of oblate JR ydenl 
19921 lLambas et alJ Il992t iTremblav fc Merrittl 119951 
Alam fc Rvdenl 12003) or prolate iTremblav fc Merritti 
1995) spheroids. iMerritt fc Tremblavl l)199fit) presented 
intriquing evidence that there exists a bimodal distri- 
bution of shapes, namely low luminosity ellipticals are 
consistent with being oblate spheroids, while bright el- 
liticals are only consistent with a triaxial intrinsic shape 
that is rounder on average than low-l uminosity ell i ticals . 
A similar conclusion was reached bv lFranx et al.l l|1991f) 
from analysis of a much smaller set of ellipticals that also 
contained kinematic information. In their study, the tri- 
axial shapes are necessary to match the small percentage 
of galaxies with rotation observed along their photomet- 
ric minor axis (see § 13.71 for additional details) and we 
note that their mean triaxiality parameter T was less 
than 0.4. We find that the mean triaxiality parameter 
for the dissipational remnants is 0.28 as opposed to 0.55 
for the dissipationless remnants. 



literature). 

The axis ratios b and c are defined as b = (A2/A1) 1 / 2 
and c = (A3/A1) 1 / 2 , where Ai > A2 > A3 are the eigen- 
values of the inertia tensor. Values of b, c, al ong with the 
"triax iality" parameter! 1 = (1 — 6 2 )/(l-c 2 ) ijFranx et all 
Il991|) are listed in Table [3] for all of the dissipationless 
and dissipational remnants. For reference, T = 1 is pro- 
late, and T = is oblate, and all values in between are 
triaxial. 

While Table |3 is useful for assessing the effects of dis- 
sipation for any particular merger, the overall trends are 
more apparent if one plots the axial ratio a against b 
as is done in Figure 1101 or a histogram of the triaxial- 
ity parameter T as in Figure ^2 Both of these figures 
show that dissipationless remnants are triaxial with a 
tendency to be more prolate. This result is also appar- 
ent in previous studies of dissipationl ess mergers ijBarnesl 
119921 lHernauisdll992t lSpringe]ll2000jK 

Figures 1 1 01 and ITTI demonstrate that dissipational rem- 
nants are also triaxial, however these remnants tend to 
be much closer to oblate. This result is consis tent wit h 
previ ous simulations that con tained 20% gas lS_prina 
2000) , or a dense stellar bulge |Hernauistlll993tlSprin 



2000; tcTonzalez-ftarcia fc Balcellsl I2005D . It is notewor- 
thy that the three dissipational systems that appear to 
be oblate isotropic rotators (d, /, and i) from their swath 
on the anisotropy diagram are actually quite triaxial 
(T = 0.34,0.41, and 0.24, respectively). However, the 
shape of merger remnants is known to be a function of 
radius ijBarne a lll)9ltlHerrIa^ and by sim- 

ply repeating our analysis on particles which are progres- 
sively more bound, we find that the central regions of the 
remnants are much more oblate than the outer regions. 
Therefore, the apparent discrepancy between the rota- 
tion of the remnant and its shape is likely a byproduct 
of our selecting the most bound half of stellar particles 
rather than a genuine sign of anomalous kinematics. In 
§ 13.91 we investigate how our measure of the triaxiality 
depends on orbit, gas fraction and progenitor mass. 

Because real ellipticals can be observed from only 
a single projection, their intrinsic shapes are difficult 



3.7. Minor Axis Rotation 

A number of elliptical galaxies are o bserved to have ro- 
tation along their minor axes (see e.g. JFranx et al.l|l991l 
and references therein). The presence of minor axis ro- 
tation is possible in a triaxial or prolate ellipsoid and 
impossible in an oblate, axisymmetric one ljBinnevlll985t 
Ide Zeeuw fc Franxlll99"ll) . The previous section demon- 
strated that all of our merger remnants show some degree 
of triaxiality and hence it is possible that any, or all, or 
our remnants could have some degree of minor axis rota- 
tion. However, because the dissipationless remnants are 
much more prolate (and definitely more triaxial) than 
dissipational ones, it is likely that dissipationless rem- 
nants will have more minor axis rotation. 

In order to determine the degree of minor axis rotation, 
we repeat the analysis outlined in S 12. 21 with the slit now 
placed along the minor axis of each remnant. Thus, for 
each remnant we have a rotation speed along the major 
(Vmaj) and minor (V^n) axes. We list the minor axis ro- 
tation for each of our merger remnants in Table 01 As an 
additional me asure of the minor axis rotation we follow 
Binncv (1985) and measure 



V m 



(V 2 ■ + V 

\ mai 1 » 



2 U/2' 
min/ 



(5) 



the minor axis rotation parameter. An additional mea- 
sure of the minor axis rotation, and one that is useful 
because it can be compared directly to observations, is 
the kinematic misalignment 'J, defined by 



tan* 



V m 



(6) 



[a and *, as well as the 
a , are listed in Table HJ for 



Both minor axis parameters, 
minor axis rotation speed V m i 
each merger remnant. 

Surprisingly, Table0]shows that the dissipationless and 
dissipational remnants have nearly identical minor axis 
rotation speeds. However, because the major axis rota- 
tion is much larger for the dissipational remnants, both 
fi and * are much smaller for the dissipational merger 
remnants. 
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Fig. 10. — The shape diagram for merger remnants, b and c are defined as (A2/A1) 1 / 2 , and (A3/A1) 1 / 2 , respectively. Ai, A2, and A3 are 
the eigenvalues of the inertia tensor in ascending order. Only the most bound half of particles are used when computing the inertia tensor. 
The solid line which denotes c = b would be a prolate spheroid with T = 1, where T = (1 — b 2 )/(l — c 2 ) is the triaxiality parameter. The 
right vertical axis, at b = 1 would be an oblate spheroid with T = 0. The dashed line shows T = 2/3, and the dotted line shows T = 1/3. 




TABLE 4 
Minor Axis Kinematics 



0.4 0.6 
Triaxiality, T 

Fig. 11. — Histogram of remnant "triaxi ality" parameters T , 
where T = (1 - 6 2 )/(l - c 2 ), as defined by IFranx etaT] 119911) . 
Oblate galaxies have T = 0, prolate galaxies have T = 1, and all 
values in between, as the remnants presented here, are triaxial. 




Note. — Minor axis rotation properties of the simulated 
merger remnants. Vmin is the minor axis rotation speed, in 
km s — 1 , fi is the minor axis rotation parameter, defined by 
Equation^] and ^ is the kinematic misalignment defined by 
Eauationlbl 



Figure 1121 presents another way to visualize the rela- 
tion between minor and major axis rotation by plotting 
the minor axis rotation parameter fi against the "total" 
rotation speed Vtot = (Knaj + Knin) 1 ^ 2 - AU merger rem- 
nants inhabit a clear band in the /i-Vtot plane. However, 
the dissipational and dissipationless remnants are con- 
centrated in very different regions. Fast rotators, those 
with Vtot > 50 km s , have very little minor axis rota- 
tion and are produced solely by dissipational mergers. 
Slow rotators, those with total rotation speeds below 
50 km s have a wide variety of minor axis rotation pa- 
rameters. Also, there is a clear preference for dissipa- 



tionless merger remnants to have significant minor axis 
rotation relative to the total rotation. In fact, the ma- 
jority of the dissipationless merger remnants have more 
minor axis rotation than major axis rotation. 

To confirm the large difference between the minor axis 
kinematics of dissipationless versus dissipational merger 
remnants, we show the histograms of the kinematic mis- 
alignment parameter \l/ from all projections in Figurc lTSl 
Again, a clear dichotomy in \l/ is shown. Dissipation- 
less remnants have significant minor axis rotation and 
thus large kine matic misalignment s. Overplotted are 
data compiled bv lFranx et alJ l)1991l) shown with a thick 
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Fig. 12. — The total rotation velocity plotted against the minor axis rotation parameter, fi, as defined by Equation^ The dashed vertical 
line denotes equivalent rotation along the major V ma j and minor V m i n axes. 
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Kinematic Misalignment, ^ 



Fig. 13. — Histogram of the kinematic misalignment angle \P. As 
with the previous figures, the filled histogram is for dissipational 
remnants and the open histogram is for dissipationle ss remnants. 
Overp lotted is data from 44 ellipticals compiled by IFranx et al.l 
1T9OT) . 



black histogram. While it is clear that there are elliptical 
galaxies with observed minor axis rotation, the majority 
(77%) have <3> < 20°. There seems to be a large peak of 
galaxies with little or no observed minor axis rotation, 
'A 7 ~ 0, and another, much smaller, peak at 'J ~ 90, al- 
though the latter contains only four galaxies, and thus 
its significance is difficult to assess. In any case, the ob- 
served distribution of kinematic misalignments is much 
better represented by the remnants produced in dissipa- 
tional mergers. 

3.8. Isophotal Shapes 



The isophotes of elliptical galaxies are not perfect el- 
lipses. The distortions from a perfect ellipse can be quan- 
tified by expanding the residuals from the fitted ellipse 
in a Fourier series, 

A?' = [akCos(k(j>) + bksin(k(j>)] , (7) 
k 

where Ar are the deviations from a perfect ellipse 
as a function of angle (p. The coefficient 04 of the 
cos (40) term measures deviations symmetric to the ma- 
jor axis. Positive values of 0,4 indicate disky is ophotes 
and negative values ind icate boxy isophotes l)Benderl 
rmsaiBender et aUllflS^ . 

As mentioned in the introduction, elliptical galaxies 
appear to com e in two types, either disky, or boxy This 
dichotomy led Kormendv & Bender (1996) to suggest a 
modified Hubble classification system in which ellipticals 
are delineated by their isophotal shapes as opposed to 
the standard ellipticity. One line of speculation suggests 
that disky isophotes are an indication of increased dissi- 
pation, and thus the dichotomy in isophotal shapes rep- 
resents varying amo unts of dissipation during a merger 
l)Bender et al.lTl992T) . However, numerical simulations 
have shown that remnants formed from dissipationless 
merging can appear both boxy and disky, depending 
on th e viewing angle and isophotal radius (Hcvl et al. 
1994). It has also been shown that there is a w eak trend 
for dissipational remnants to be more disky ijSpringel 
120001 iBekki k. Shiovalll997j) probably owing to the pro- 
ducti on of central density cusps that destabilize box or- 
bits llGerhard fc Binnevlll985l iMerritt fc FridmaiJll996l 
iBarnes fc HernauistJll996|) . 

We measure the deviations from perfect ellipses for 
each of our remnant images and list the average for each 
orientation in Table El From these averages we confirm 
that the dissipational remnants are more disky, in gen- 
eral, than their dissipationless counterparts. We do cau- 
tion, however, that the scatter within any particular rem- 
nant is large, reflecting that remnants can appear disky 
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TABLE 5 
Isophotal Shapes 



ID 


Dissipationlcss 
100(24 /a 


40% Gas 
100a4 /a 


h 


-1.5 


0.9 




0.3 


1.2 




-1.8 


1.6 


d 


-1.3 


1.8 


e 


-1.0 


0.4 


f 


1.0 


0.9 


g 


2.5 


0.9 


i 


0.1 


1.9 


j 


0.5 


0.1 


k 


-2.5 


0.8 


1 


2.7 


-0.3 


m 


0.4 


-0.8 


n 


0.0 


1.6 


o 


-0.9 


1.7 


P 


-0.2 


1.0 



Note. — Average isophotal shape 
parameter a^/a multiplied by one 
hundred for each simulated merger 
remnant. Positive values of a^/a de- 
note disky isophotcs, while negative 
values denote boxy isophotcs. 
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40% gas 




-6 -4 -2 2 4 6 
100 a4/a 

Fig. 14. — Histogram of 100a4/a, which measures the isopho- 
tal deviations from a perfect ellipse; positive values indicate disky 
isophotes and negative values indicate boxy isophotes. As with 
the previous figures the filled histogram is for dissipational merger 
remnants and the open histogram is for the dissipationless merger 
remnants. O verplo tted is data from 59 ellipticals compiled by 
IBender et all Il99^ . 



or boxy depending on the viewing angle. 

In Figure El we show a histogram of the coefficient 
d4 (normalized by the semi-major axis a and multiplied 
by an arbitrary factor of one hundred). Also included 
in this figure is da ta from 59 ellipticals compiled by 
IBender et alJ l|1992t K In general, there is good agreement 
between the simulated merger remnants and the data. 
While most isophotes having relatively small deviations 
from perfect ellipses, dissipational remnant tend to be 
slightly disky and dissipationlcss remnants slightly boxy. 
Observationally, there are an equal number of disky and 



boxy remnants. 

We note that there exists a tail of disky dissipational 
remnants and boxy dissipationless remnants not present 
in the data. This discrepancy may be a result of resolu- 
tion, as degrading the resolution of our image, or smooth- 
ing it, both reduce the measured isophotal deviation. We 
also note that many of the disky dissipational remnants 
are the fast rotating systems viewed edge-on. 

As evide nce for the dichotomy betwee n disky and boxy 
ellipticals, iKormendv fc Bender! l)1996j) presented corre- 
lations between the a 4 parameter and (V/a)* and minor 
axis rotation parameter \i. These correlations are pre- 
sented in Figure El Neither series of remnants presents 
a clear dichotomy between disky/boxy and rotation, as is 
seen in the data. The dissipationless remnants, in partic- 
ular, appear to be an especially poor match to the data 
as there are a significant number of slow, (V/a)* <C 1, 
and minor axis, fi ~ 1, rotators with disky isophotes. El- 
liptical galaxies with these rotational properties are all 
boxy, in disagreement with the simulations. 

Overall, the dissipational remnants are a much better 
match to the general trends present between the rotation 
and isophotal shape of elliptical galaxies. Most of these 
remnants have very little minor axis rotation, fj, < 1, 
(V/a)* « 1, and have isophotes that are both disky and 
boxy. There are a large number of observational analogs 
to these systems. However, Figure El makes it evident 
that we are still unable to produce slow and minor axis 
rotators that are uniformly boxy. This will be discussed 
further in § 0] 

3.9. Dependencies 

The previous sections demonstrate that there is a sig- 
nificant difference between the size, shape and rotational 
properties of merger remnants when the initial disk con- 
tains 40% gas compared to when the entire disk is dis- 
sipationless. However, it is unclear whether these re- 
sults are sensitive to the various assumptions implicit to 
our merger simulations. While a complete exploration of 
the large parameter space that describes the disk galaxy 
models and their encounters is beyond the scope of this 
work, we can provide some indications of the robustness 
of our results by systematically varying several param- 
eters which we suspect may affect the properties of the 
merger remnants. Specifically, the following three sec- 
tions consider variations in the progenitor disk gas frac- 
tion, the orbital angular momentum, and the size (mass) 
of the progenitor disk galaxies. 

For simplicity we select the e, h, and k disk orienta- 
tions from Table^to resimulate throughout this section. 
These three are selected to span a range of remnant prop- 
erties. 

3.9.1. Dependence on Gas Fraction 

To begin, we systematically change the disk gas frac- 
tion from zero, i.e., the dissipationless case, to a pure gas 
disk. For reference, our fiducial dissipational simulation 
contained a disk with 40% gas. In FigurelTfilwe show the 
size, as measured by the semi-major axis, ellipticity, and 
(V/a)*, all averaged over all projections for each merger 
remnant. In addition, we show the triaxiality parameter 
T, and the anisotropy S for each remnant. 

As might be expected from the results of § EI there is a 
systematic trend for remnants produced by the collision 



0.5 
0.0 

> -0.5 
60 

o 

J -1.0 

-1.5k 
1.0 - 



sr; 0.8 
+ 0-6 k 



> s 0.4 



1 0.2 
> 

o.o k_ 

-2 




Kinematic Structure of Merger Remnants 

0.5 



17 



• 


I 


• 

• 

• 


•I 

1 


• • •• 
• 


I 

I 

■ - 


• 
• 

• • 




• • 




•• 
m 


■ UM ** 

i * 



1 
100 a 4 /a 




1 
100 a 4 /a 



Fig. 15. — A density map that displays correlations between isophotal shape and rotation. 100a4/a, which measures the isophotal 
deviations from a perfect ellipse; positive values indicate disky isophotes and negative values indicate boxy isophotes. The top panel shows 
the rotation parameter (V/cr)* and the l ower panel shows the minor axis rotation parameter of Eq. [5] Overplotted is the identical data as 
presented in Kormcndy & Bender (1996). 



of spiral galaxies with an increasing fraction of gas to be 
smaller, rounder, and to have more rotation. However, 
some of the remnant properties are more sensitive to the 
gas fraction than others. For example, the size of rem- 
nants shrinks and (V/cr)* steadily rises with increasing 
gas fraction. But, the projected ellipticity becomes only 
slightly rounder, and at a level that is comparable to the 
variations expected from projection effects alone (shown 
by the error bars on each point). For each individual 
orientation the shape, as measured by the triaxiality pa- 
rameter T, appears to level off above a gas fraction of 
20%, even though the individual orientations level off to 
different shapes. 

The dependencies displayed by trends in (V/cr)*, T, 
and S demonstrate that inherent differences exist purely 
because of the initial disk orientations. For example, 
the e merger tends to be much slower rotating, much 
more triaxial, and contain much less anisotropy, no mat- 
ter what the gas fraction is, even though the trends with 
gas fraction are similar to those found in h and k mergers. 

From this series of 21 simulations (3 orientations, each 
simulated with 7 different gas fractions) we conclude that 
the size and rotation are strongly affected by the initial 
disk gas fraction, while the shape is only moderately af- 
fected. The anisotropy appears more strongly correlated 
with the orbital configuration than with the gas fraction. 

These results also directly relate to those of 
iRobertson et al.l l)2006a|) who find that remnants pro- 
duced in major mergers can reproduce the observed fun- 
damental plane of ellipticals when the progenitor disks 
contain gas fractions greater than 30%. Thus, it seems 
that gas-rich major mergers can help explain both the 
kinematics of elliptical galaxies and the tilt of the funda- 
mental plane. 

3.9.2. Dependence on Orbital Angular Momentum 

Next, we note that merger orbits as determined by 
cosmological simulations span a range of impact param- 



eters (Khochfar & Burkert 200 lj. In the fiducial series 
of mergers, analyzed in § [21 we assumed each merger 
orbit was parabolic, and had a pericentric distance of 
7.1 kpc, generating an orbit with a low impact parame- 
ter. In this series of runs, we resimulate the dissipational 
mergers (e, h, and k) with a range of pericenter separa- 
tions between 3.6 and 57 kpc. The initial starting sepa- 
ration for all mergers was 143 kpc. We note that all of 
these orbits contain an identical amount of energy (zero) 
yet have progressively more orbital angular momentum 

[L ~ -Rp eri ) ■ 

Figure El is equivalent to Figure El except in the 
present figure the horizontal axis measures the pericen- 
tric distance. The size of each merger remnant is a mod- 
erate function of i? pC ri, with wide orbits producing larger 
remnants. The apparent ellipticity has no discernible 
correlation with i? pcr i- The rotation shows weak, or no, 
correlation to i? pC ri- At first glance, this lack of corre- 
lation between rotation and orbital angular momentum 
seems odd. However, two effects act to offset the transfer 
of orbital angular momentum to remnant rotation. First, 
the initial disk galaxies contain massive dark halos which 
soak up the majority of the orbital angular momentum. 
Second, the orbits with high orbital angular momentum 
require a longer time to merge, and thus the initial disks 
have consumed a larger fraction of their initial gas. In 
this sense, they are effectively lower gas fraction mergers, 
which the previous section showed are slower rotators. 

As we saw in the last section, the shapes of merger rem- 
nants depend on the original disk orientations. Figure lT7l 
demonstrates that the shape also depends on i? pcr i. Wide 
orbits tend to produce prolate (T w 1) remnants, regard- 
less of the initial disk orientation, while intermediate and 
radial orbits produce remnants that are closer to oblate. 
However, the orbits with i? pcr i < 30 kpc begin to segre- 
gate based upon the initial disk orientation. There is a 
trend, though, for very radial orbits to be more prolate 



18 



Cox, et al. 



0. 

A 

v 



V 



E- 1 



E- 1 



'o 



0.6 



2 o 

o 
< 



0.0 



X— * orbit e 
orbit h 







0.0 0.2 0.4 0.6 0.8 1.0 
Progenitor Gas Fraction 



Fig. 16. — The semi-major axis a, ellipticity e, and V ma j/(r aver- 
aged over 190 different lines of sight and the triaxiality parameter 
T, and anisotropy parameter S for our e, h (prograde-prograde), 
and k orientation merger remnants as a function of the gas fraction 
of the progenitor disk. A key is provided in the top panel which 
denotes the symbols used for the three initial disk orientations. For 
the top three panels, which present quantities averages over many 
lines of sight, one sigma error bars are included to provide a sense 
of the variation introduced by projection effects. 



than intermediate orbits. 

There is a mild tendency for very wide orbits to pro- 
duce more isotropic (<5 ~ 0) remnants. This correlation, 
however, is weak, at best, especially compared to the 
strong dependence on the initial disk orientations. 

3.9.3. Dependence on Progenitor Mass 

As a final dependency we investigate the relation 
between progenitor mass and the remnant properties. 
In order to span a wide range of remnant central 
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Fig. 17. — Same as Figure lTrTl exccDt here the remnant properties 
are plotted against the orbital pericentric distance. 



velocity dispersions, we have simulated our standard 
e, h, and k mergers with initial disk galaxies both 
smaller (V200 = 56,80,115 km s" 1 ) and larger (V200 — 
220,320,500 km s" 1 ) than our fiducial case (V200 — 
160 km s _1 ). For these mass excursions, all other pro- 
genitor disk properties are kept fixed, such as gas fraction 
(40%), spin (0.033), dark matter concentration (9), and 
specific angular momentum (0.05). The initial orbits are 
scaled such that the ratio R pC ii/Rd remains constant. 
For reference, our fiducial model (V200 = 160 km s _1 ) 
has a total virial mass of 1.4 x 10 12 Mq, and our mass 
excursion samples galaxies a factor of 20 smaller and a 
factor of 30 larger in mass. 

Figure ITH1 shows the resulting remnant properties as a 
function of remnant central velocity dispersion. In gen- 
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Fig. 18. — Same as Figure UHl exccpt here the remnant properties 
are plotted against the central velocity dispersion a of each merger 
remnant. 



eral, it appears as if the size is the strongest dependency, 
as expected since Rd oc Vion- We find a size-mass rela- 
tionship consistent with lRobertson et alJ l)2006bla(h who 
simulate binary mergers for an equivalent range of pro- 
genitor disk masses. The ellipticity and rotation do not 
appear to correlate with mass. Once again, the triaxial- 
ity and anisotropy depend on the initial disk orientations 
more than on the size of the merging systems. 

From our investigation of the remnant properties as a 
function of progenitor mass we conclude that mass is of 



little importance in fixing the rotational properties of a 
merger remnant. Instead, it is the progenitor gas fraction 
and merger orbit and orientation which are the primary 
determinants in what is left behind from the collision of 
two disk galaxies. 

3.9.4. Other Dependencies 

In the three previous subsections we have tried to de- 
termine the degree to which our results are dictated by 
our initial assumptions. We conclude that the initial gas 
fraction is of critical importance to the structure of a 
merger remnant. The orbit and orientation, as well as 
the progenitor disk mass are also factors which deter- 
mine the remnant properties. However, there are still a 
large number of parameters (bulge size, radial gas pro- 
file, halo A and c) which may in principle influence the 
remnants that have not been addressed so far. 

We have performed a single simulation where the initial 
disk contained a bulge one-third the mass of the progen- 
itor disk. The fiducial disk galaxy model did not con- 
tain a bulge. In this case the remnant showed system- 
atic, but very small differences from the no-bulge fiducial 
remnant. Specifically, because of the dense central bulge, 
the remnant is rounder and has a higher central velocity 
dispersion. The no-bulge remnant also has more rota- 
tion, but because a increases, the net effect is to lower 
V/a, and drive the remnant down and to the left in the 
anisotropy diagram. The presence of a bulge increases 
the rotation because star formatio n is suppressed dur- 
ing the earl y stages of the merger (Mih os fc HernauisH 
Il994d Il996j) and hence the effective gas fraction at the 
final merger is larger than in the no-bulge case. We ex- 
pect that any process, or assumption, that results in a 
high gas fraction during the final merger will increase the 
rotation of the remnant. Possible mechanisms are the in- 
clusion of a bulge, less efficient quiescent star formation, 
or an extended gas distribution. 

We finally note that the presence of a black hole has 
very little effect to the results presented here. The black 
hole feedback induced "blowout" phase is necessary to 
termin ate star formation an d thereby produce red rem- 
nants l)Springel et al.ll2005a|) . but this process occurs af- 
ter the merger has taken place and subsequent to the 
production of most of the remnant stars. In fact, we ran 
one case that did not include black hole feedback, and 
thus did not produce a galactic wind, however the rem- 
nant stellar mass increased by less than 5%. While these 
additional stars add to the "young" component shown in 
Figure the overall changes to the remnant kinematics 
are much less than projection effects or merger orien- 
tations. Thus we conclude that the black hole may be 
very important for the long-term color evolution, but has 
little importance in determining the bulk kinematics of 
the remnant because the processes that fix the dynamics 
have finished before the blowout. 

4. DISCUSSION 

The previous section presented an analysis of two se- 
ries of fifteen merger simulations, one series consisted 
of dissipationless mergers, while the second series con- 
tained disks composed of 40% gas, but identical in every 
other respect. In addition, we investigate the dependence 
of our results to variations in the original gas fraction, 
orbital angular momentum, and progenitor size/mass. 
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These tests show that the size, shape and rotation of any 
merger remnant is a strong function of the initial disk gas 
fraction and the orientations of the disks. In what fol- 
lows we will attempt to place these results into a broader 
context and relate them to our current understanding of 
the formation and evolution of spheroidal galaxies. 

4.1. Forming Low-luminosity Elliptical Galaxies and 

Bulges 

As outlined in the Introduction, low-luminosity ellip- 
ticals and bulges fall into the class of galactic spheroids 
that have disky isophotes, are X-ray faint, have power- 
law surface brightness profiles and show significant ro- 
tation. The rotation axis is typically aligned with their 
photometric minor axis and they are consistent with be- 
ing oblate isotropic rotators. Thus, asking "What mech- 
anism forms low-luminosity ellipticals and bulges?" may 
be equivalent to the question "What produces oblate 
isotropic rotators?" 

Our results suggest that the merger of equal mass gas- 
rich disk galaxies may produce remnants that are oblate 
isotropic rotators. Qualitatively, over 25% (4/15) of the 
40% gas major mergers yielded an oblate isotropic rota- 
tor. However, it is unclear if major mergers can produce 
these objects in the correct number. As indicated by 
FigureEl t ne dissipational mergers produce slightly fewer 
fast rotating, (V/a)* > 1, ellipticals than observed. Does 
this dearth of fast rotating remnants pose a problem to 
the "merger hypothesis" ? 

The short answer to this question is no. The sce- 
nario outlined in the current paper, major mergers be- 
tween gas-rich disk galaxies, is likely only one of a num- 
ber of possible mechanisms to produce oblate isotropic 
rotators. Another po ssible mechanism, proposed by 
iNaab fc Burkerd (J2003) , is the dissip ationless merger 
of une qual mass disk galax i es. Both [ Naab fc BurkerU 
( 2003) and iBournaud et all <|2004 1200.51) show that un- 
equal mass mergers, specifially those with mass ratios 1:4 
or higher, produce oblate remnants. 

While we consider there to be strong evidence that 
all mergers should include progenitor galaxies with some 
amount of gas, we suspect that the inclusion of gas 
will only increase the amount of rotation in the rem- 
nants of minor mergers. Finally, we note that Figure ITH1 
shows that the frequency of oblate isotropic rotators will 
increase when initial disks of higher gas fraction are 
merged. In short, it does not apear that there are any 
problems in producing oblate isotropic spheroids. 

It is also noteworthy that not all of the gas is con- 
sumed by the merger event. In general, ~ 10 — 20% 
of the original gas mass remains in the remnant, typi- 
cally in a hot phase that is spread throughout the dark 
halo ijCox et alj[2004l [20061 The cooling time of this 
gas is several Gyr, and thus, depending on the angular 
momentum of this gaseous material, and the amount of 
pristine inter-galactic material that is newly accreted, it 
is possible a new gas disk will from around, or within, 
the existing spheroidal stellar remnant. In this sense the 
remnants produced here may be the precursor to bulges 
in present day spirals or at least a hidden disk compo- 
nent. 

4.2. Forming Large Elliptical Galaxies 



Luminous elliptical galaxies have boxy isophotes, con- 
tain halos of hot gas, have surface-brightness profiles with 
cores and demonstrate little or no rotation. While form- 
ing low-luminosity ellipticals was equivalent to forming 
oblate isotropic rotators, an analogous equation does not 
exist for luminous ellipticals. We do note that the pri- 
mary conclusion of our work is that dissipational merger 
remnants are a much better match to the entire class 
of elliptical galaxies. The evidence for this comes from 
the distribution of ellipticities (Figure |3J) and (V/a)* 
(Figure EJ), and the amount of minor axis rotation (Fig- 
ure El j as measured by 'J, the kinematic misalignment 
as well as the correlations between rotation and isopho- 
tal shape (Figure IT31) . However, as the isophotal analy- 
sis demonstrates, it is unclear that any of our remnants, 
dissipational or dissipationless, are really slow rotators 
and uniformly boxy as it appears all luminous ellipticals 
are. This difficulty with the model brings into question 
whether luminous ellipticals can be produced within the 
merger hypothesis. 

While a single gas-rich merger appears to well re- 
produce the properties of low- luminosity ellipticals, 
it is possible that luminous ellipticals require a more 
complicated merger history. In fact, within the hierar- 
chical build-up of structure predicted by Lambda cold 
dark matter, we expect mergers to occur frequently, 
especially minor mergers, and there is no reason to 
believe that galaxies participate in only one major 
merger. Recent observations suggest that, on average, 
somewhere between one-half and al l elliptical galaxies 
undergo a spher oid-spheroid merger l|van Dokkuml l2005: 
iBell et al J 1200^1 below redshift one. In the context of 
forming bright cluster galaxies, multiple dry mergers 
has l ong been a we ll-motivated mechanism (jMerritti 
119851: iDubinski [199^1 . As a brief test of the effects of 
continued merging, we performed one re-merger, where 
two merger remnants were merged in a manner similar 
to our fiducial series. In this one case, the remnant rota- 
tion and ellipticity decreased slightly, while the central 
velocity dispersion increased slightly and the isophotes 
became increasingly boxy. Although it is difficult to say 
anything concrete from a single simulation, we speculate 
that spheroid-spheroid merging will move remnants 
vertically downwards in the anisotropy diagram to the 
region occupied by luminous slowly-rotating ellipticals 
while also making the system uniformly boxy. A similar 
result has been found for m ergers of dissipationless 
remnants bv INaab et alJ ||200(tD and we note that other 
work has also shown that multiple m erger remnants 
can remain on the fundamental plane lC^rjeialQ_£LaI 
1995[ iNmoti et al.ll2008t lOonzalez-Garcia fc /jU^AIbada 
20031 "iBovlan-Kolchin et alJ l2006t iRobertson et al. 
2006al) and may have i? 1 / 4 surface brightness profiles 
ijBournaud et alj|2005jl . However, it is clear that much 
more work needs to be performed in order to understand 
the relationship between different merger histories and 
the kinematics of the remnant. 

4.3. Future Considerations 

As a final point of discussion we note that a consid- 
erable amount of work still needs to be performed in 
order to prove that the simulated merger remnants are 
bona fide elliptical galaxies. The present study has shown 
that dissipational merger remnants appear to be a good 



Kinematic Structure of Merger Remnants 



21 



match (at least much better than dissipationless rem- 
nants) to the cllipticities, (V/cr)*, and kinematic mis- 
alignments observed in elliptical galaxies. However, 
we have not addressed the surface brightness profiles, 
metallicities and abundance ratios, and kinematic sub- 
systems of our merger remnants. Work toward this end 
is currently underway and is necessary as many of these 
features are observed to be correlated with the rotational 
properties suggesting that all of these properties are a re- 
sult of a common formation mechanism. 

5. CONCLUSIONS 

In this paper we have analyzed a series of numer- 
ical simulations to investigate the kinematic structure 
of galaxies formed from the collision of equal-mass disk 
galaxies. In particular, we determine that remnants 
produced by the collision of gas-rich disk galaxies are 
smaller, rounder, have higher central velocity dispersions 
(on average) and larger rotation speeds than the rem- 
nants of dissipationless disk galaxy mergers. 

The larger rotation present in gas-rich merger rem- 
nants owes its origin to dissipation. Stars formed during 
the merger rotate faster than stars present before the 
merger began. Dissipation and star formation also pro- 
duce remnants that are closer to oblate and are uniformly 
more isotropic than dissipationless remnants. Slightly 
more than one-quarter of the dissipational merger rem- 
nants are consistent with being oblate isotropic rotators. 

Compared to observed ellipticals, the remnants formed 
from dissipational mergers are a significantly better 
match to the distribution of ellipticities, (V/cr)*, and 
kinematic misalignment angles $ of elliptical galaxies 
than dissipationless merger remnants. In particular, dis- 
sipationless remnants demonstrate significant minor axis 
rotation and appear to be flattened much more than ob- 
served ellipticals. 

We also calculate the isophotal shapes of the simu- 
lated merger remnants. Dissipational remnants tend to 
be disky while dissipationless remnants tend to be boxy. 
Observed ellipticals are evenly distributed between disky 
and boxy. While both remnants appear to be a sufficient 
match to the observed distribution, when comparing the 
correlation between disky/boxy and rotation, the dissi- 
pationless merger produce a significant number of slowly 
rotating disky remnants where there are no observed 
analogs. Both dissipational and dissipationless mergers 
produce remnants that are both disky and boxy, and thus 
these mechanisms have difficulty reproducing the lumi- 
nous, slowly rotating ellipticals that are observed to be 
uniformly boxy. 

In general, our results suggest that dissipationless disk 
galaxy mergers cannot be the dominant mechanism to 



form elliptical galaxies. Dissipational mergers, on the 
other hand, appear to be a viable mechanism to pro- 
duce elliptical galaxies, specifically oblate isotropic ro- 
tators (i.e., low-luminosity ellipticals), and thus our re- 
sults lend support to the "merger hypothesis" provided 
that the progenitor disk galaxies are gas-rich. As men- 
tioned in § ^ additional evidence for the dissipative ori- 
gin of ellipticals comes from their high phase space den- 
sity compared to spiral galaxies. Gas dissipation pro- 
vides a natural mechanism to increase the phase space 
density during the merger and also appears to be nec- 
essary fo r reproducing the scali ng relation of elliptical 
galaxies l|Robertson et alJl2006al. 

Our modeling suggests several avenues for further 
testing this hypothesis. Mergers between gas-rich spi- 
rals will imprint subtle features into the remnants. 
The central st arburst will modify the in ner profiles of 
the remnants (Mi hos fe Hernauistl Il994a|) . perhaps ac- 
counting central light excesses s een in merging systems 
(|Rothberg fc JosenrJ 120041 120061. In principle, this can 
be tested by comparing predictions for metallicity and 
color gradients with observations (Mih os fc Hernauistl 
1994b). Dissipational merger may also provide a natural 
mechanism to produce kinematic subsystems in ellipti- 
cal ga laxies l|Hernauist fc Barneslll99lHBendo fc Barnes! 
2000). The shells, ripples, loops and other fine struct ures 
seen around many relaxed ellipticals ( Schwciz erll998fl are 
a n atural consequence of merg ers involving disk galax- 
ies l)Hernauist fc Snergellll992fh but that do not form in 
major mergers between hot stellar systems. Determin- 
ing the ubiquity of fine structure in red galaxies would 
further constrain the importance of disk mergers to the 
formation of ellipticals. 

Placed within the "cosmic cycle" of galaxy formation, 
we can now argue that gas-rich major mergers trigger 
quasars and starbursts, fuel the growth of supcrmas- 
sive black holes, and produce remnant galaxies which 
have the colors, scaling relations, and kinematics akin to 
present day ellipticals. 



We thank Marijn Franx for the motivation that ini- 
tiated this work and additional comments that signifi- 
cantly improved this paper. We also thank Ralf Bender 
for providing the data in Figure 1151 This work was sup- 
ported in part by NSF grants ACI 96-19019, AST 00- 
71019, AST 02-06299, and AST 03-07690, and NASA 
ATP grants NAG5-12140, NAG5-13292, and NAG5- 
13381. The simulations were performed at the Cen- 
ter for Parallel Astrophysical Computing at Harvard- 
Smithsonian Center for Astrophysics. 



REFERENCES 



Alam, S. M. K. & Ryden, B. S. 2002, ApJ, 570, 610 

Allgood, B., Flores, R. A., Primack, J. R., Kravtsov, A. V., 

Wechsler, R. H., Faltenbacher, A., & Bullock, J. S. 2006, 

MNRAS, 367, 1781 
Barnes, J. E. 1988, ApJ, 331, 699 
— . 1992, ApJ, 393, 484 

Barnes, J. E. & Hernquist, L. 1996, ApJ, 471, 115 
Barnes, J. E. & Hernquist, L. E. 1991, ApJ, 370, L65 
Bckki, K. & Shioya, Y. 1997, ApJ, 478, L17+ 



Bell, E. F., Naab, T., Mcintosh, D. H., Somerville, R. S., Caldwell, 
J. A. R., Barden, M., Wolf, O, Rix, H.-W., Beckwith, S. V. W., 
Borch, A., Haeussler, B., Heymans, O, Jahnke, K., Jogee, S., 
Meisenheimer, K., Peng, C. Y., Sanchez, S. F., & Wisotzki, L. 
2005, ArXiv Astrophysics e-prints 

Bender, R. 1988, A&A, 193, L7 

Bender, R., Burstein, D., & Faber, S. M. 1992, ApJ, 399, 462 
Bender, R., Doebereiner, S., & Moellenhoff, C. 1988, A&AS, 74, 
385 

Bender, R. & Nieto, J.-L. 1990, A&A, 239, 97 



22 



Cox, et al. 



Bender, R., Surma, P., Doebereiner, S., MoellenhofT, C, & 

Madcjsky, R. 1989, A&A, 217, 35 
Bendo, G. J. & Barnes, J. E. 2000, MNRAS, 316, 315 
Binney, J. 1978, MNRAS, 183, 501 
— . 1985, MNRAS, 212, 767 

Binney, J. & Tremaine, S. 1987, Galactic dynamics (Princeton, NJ, 

Princeton University Press, 1987, 747 p.) 
Bournaud, F., Combes, F., & Jog, C. J. 2004, A&A, 418, L27 
Bournaud, F., Jog, C. J., & Combes, F. 2005, A&A, 437, 69 
Boylan-Kolchin, M., Ma, C.-P., & Quatacrt, E. 2005, MNRAS, 362, 

184 

— . 2006, astro-ph/0601400 

Burkert, A. & Naab, T. 2005, MNRAS, 363, 597 
Capelato, H. V., de Carvalho, R. R., & Carlberg, R. G. 1995, ApJ, 
451, 525 

Cappellari, M., Bacon, R., Bureau, M., Damcn, M. C, Davics, 
R. L., de Zeeuw, P. T., Emscllem, E., Falcon-Barroso, J., 
Krajnovic, D., Kuntschner, H., McDermid, R. M., Peletier, R. F., 
Sarzi, M., van den Bosch, R. C. E., & van dc Vcn, G. 2006, 
MNRAS, 366, 1126 

Carlberg, R. G. 1986, ApJ, 310, 593 

Cox, T. J., Di Matteo, T., Hernquist, L., Hopkins, P. F., Robertson, 
B., & Springel, V. 2006, ApJ accepted (astro-ph/0504156) 

Cox, T. J., Jonsson, P., Primack, J., & Somerville, R. S. 2005, 
MNRAS, submitted (astro-ph/0503201) 

Cox, T. J., Primack, J., Jonsson, P., & Somerville, R. S. 2004, ApJ, 
607, L87 

Cretton, N., Naab, T., Rix, H., & Burkert, A. 2001, ApJ, 554, 291 

Davies, R. L. & Birkinshaw, M. 1988, ApJS, 68, 409 

Davies, R. L., Efstathiou, G., Fall, S. M., Illingworth, G., & 

Schechter, P. L. 1983, ApJ, 266, 41 
de Vaucouleurs, G. 1948, Annales d'Astrophysique, 11, 247 
de Zeeuw, P. T., Bureau, M., Emsellem, E., Bacon, R., Marcella 

Carollo, C, Copin, Y., Davics, R. L., Kuntschner, H., Miller, 

B. W., Monnet, G., Peletier, R. F., & Verolme, E. K. 2002, 

MNRAS, 329, 513 
de Zeeuw, T. & Franx, M. 1991, ARA&A, 29, 239 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 

604 

Dubinski, J. 1998, ApJ, 502, 141 

Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, N. A., 

& Adelberger, K. L. 2006, ApJ accepted (astro-ph/0604041) 
Faber, S. M., Tremaine, S., Ajhar, E. A., Byun, Y., Dressier, A., 

Gcbhardt, K., Grillmair, C, Kormcndy, J., Lauer, T. R., & 

Richstone, D. 1997, AJ, 114, 1771 
Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., 

Lynden-Bell, D., & Terlevich, R. J. 1989, ApJS, 69, 763 
Ferrarese, L. & Merritt, D. 2000, ApJ, 539, L9 
Franx, M., Illingworth, G., & de Zeeuw, T. 1991, ApJ, 383, 112 
Franx, M., Illingworth, G., & Heckman, T. 1989, ApJ, 344, 613 
Gcbhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., 

Filippcnko, A. V., Green, R., Grillmair, C, Ho, L. C, Kormcndy, 

J., Lauer, T. R., Magorrian, J., Pinkney, J., Richstone, D., & 

Tremaine, S. 2000, ApJ, 539, L13 
Gerhard, O. E. & Binney, J. 1985, MNRAS, 216, 467 
Gonzalez-Garcfa, A. C. & Balcells, M. 2005, MNRAS, 357, 753 
Gonzalez-Garci'a, A. C. & van Albada, T. S. 2003, MNRAS, 342, 

L36 

— . 2005, MNRAS, 361, 1043 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, 

F. K., Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759 
Hernquist, L. 1990, ApJ, 356, 359 
— . 1992, ApJ, 400, 460 
— . 1993, ApJ, 409, 548 

Hernquist, L. & Barnes, J. E. 1991, Nature, 354, 210 
Hernquist, L. & Spergel, D. N. 1992, ApJ, 399, L117 
Hernquist, L., Spergel, D. N., & Heyl, J. S. 1993, ApJ, 416, 415 
Heyl, J. S., Hernquist, L., & Spergel, D. N. 1994, ApJ, 427, 165 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Martini, 

P., Robertson, B., & Springel, V. 2005a, ApJ, 630, 705 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B., & Springel, V. 2005b, ApJ, 630, 716 



— . 2005c, ApJ, 632, 81 
— . 2006a, ApJS, 163, 1 

Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., Di Matteo, 

T., & Springel, V. 2006b, ApJ, 639, 700 
Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Springel, 

V. 2006c, ApJS, 163, 50 
Hopkins, P. F., Hernquist, L., Martini, P., Cox, T. J., Robertson, 

B., Di Matteo, T., & Springel, V. 2005d, ApJ, 625, L71 
Kazantzidis, S., Mayer, L., Colpi, M., Madau, P., Dcbattista, V. P., 

Wadsley, J., Stadel, J., Quinn, T., & Moore, B. 2005, ApJ, 623, 

L67 

Kennicutt, R. C. 1998, ApJ, 498, 541 
Khochfar, S. & Burkert, A. 2001, ApJ, 561, 517 
Kormendy, J. & Bender, R. 1996, ApJ, 464, L119+ 
Lake, G. 1989, AJ, 97, 1312 

Lambas, D. G., Maddox, S. J., & Lovcday, J. 1992, MNRAS, 258, 
404 

Lauer, T. R., Faber, S. M., Gcbhardt, K., Richstone, D., Tremaine, 
S., Ajhar, E. A., Aller, M. C, Bender, R., Dressier, A., 
Filippcnko, A. V., Green, R., Grillmair, C. J., Ho, L. C, 
Kormendy, J., Magorrian, J., Pinkney, J., & Siopis, C. 2005, 
AJ, 129, 2138 

Lidz, A., Hopkins, P. F., Cox, T. J., Hernquist, L., & Robertson, 

B. 2006, ApJ, 641, 41 
Merritt, D. 1985, ApJ, 289, 18 
Merritt, D. & Fridman, T. 1996, ApJ, 460, 136 
Merritt, D. & Tremblay, B. 1996, AJ, 111, 2243 
Mihos, J. C. & Hernquist, L. 1994a, ApJL, 437, L47 
— . 1994b, ApJ, 427, 112 
— . 1994c, ApJ, 431, L9 
— . 1996, ApJ, 464, 641 

Milosavljevic, M. & Merritt, D. 2001, ApJ, 563, 34 
Naab, T. & Burkert, A. 2003, ApJ, 597, 893 
Naab, T., Khochfar, S., & Burkert, A. 2006, ApJ, 636, L81 
Negroponte, J. & White, S. D. M. 1983, MNRAS, 205, 1009 
Nipoti, C, Londrillo, P., & Ciotti, L. 2003, MNRAS, 342, 501 
Ostriker, J. P. 1980, Comments on Astrophysics, 8, 177 
Rix, H., Carollo, C. M., & Freeman, K. 1999, ApJ, 513, L25 
Robertson, B., Cox, T. J., Hernquist, L., Franx, M., Hopkins, P. F., 

Martini, P., & Springel, V. 2006a, ApJ, 641, 21 
Robertson, B., Hernquist, L., Bullock, J. S., Cox, T. J., Di Matteo, 

T., Springel, V., & Yoshida, N. 2005, arXiv:astro-ph/0503369 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 

P. F., Martini, P., & Springel, V. 2006b, ApJ, 641, 90 
Rothberg, B. & Joseph, R. D. 2004, AJ, 128, 2098 
— . 2006, AJ, 131, 185 
Ryden, B. 1992, ApJ, 396, 445 
Schmidt, M. 1959, ApJ, 129, 243 

Schweizer, F. 1998, in Saas-Fee Advanced Course 26: Galaxies: 

Interactions and Induced Star Formation, ed. D. Friedli, 

L. Martinet, & D. Pfenniger 
Sersic, J. L. 1968, Atlas de galaxias australes (Cordoba, Argentina: 

Observatorio Astronomico, 1968) 
Shen, S., Mo, H. J., White, S. D. M., Blanton, M. R., Kauffmann, 

G., Voges, W., Brinkmann, J., & Csabai, I. 2003, MNRAS, 343, 

978 

Springel, V. 2000, MNRAS, 312, 859 
— . 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005a, ApJ, 620, L79 
— . 2005b, MNRAS, 361, 776 

Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 
— . 2003, MNRAS, 339, 289 

Toomrc, A. 1977, in Evolution of Galaxies and Stellar Populations, 
p.401 

Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Tremblay, B. & Merritt, D. 1995, AJ, 110, 1039 

Trujillo, I., Aguerri, J. A. L., Cepa, J., & Gutierrez, C. M. 2001, 

MNRAS, 321, 269 
van Dokkum, P. G. 2005, AJ, 130, 2647 



